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Abstract 

A  Swinging  Atwood  Machine  (SAM)  is  built  and  some  experimental  results  concerning 
its  dynamic  behaviour  are  presented.  Experiments  clearly  show  that  pulleys  play  a  role 
in  the  motion  of  the  pendulum,  since  they  can  rotate  and  have  non-negligible  radii  and 
masses.  Equations  of  motion  must  therefore  take  into  account  the  inertial  momentum 
of  the  pulleys,  as  well  as  the  winding  of  the  rope  around  them.  Their  influence  is  com¬ 
pared  to  previous  studies.  A  preliminary  discussion  of  the  role  of  dissipation  is  included. 
The  theoretical  behaviour  of  the  system  with  pulleys  is  illustrated  numerically,  and  the 
relevance  of  different  parameters  is  highlighted.  Finally,  the  integrability  of  the  dynamic 
system  is  studied,  the  main  result  being  that  the  Machine  with  pulleys  is  non-integrable. 
The  status  of  the  results  on  integrability  of  the  pulley-less  Machine  is  also  recalled. 

Key  words:  Swinging  Atwood’s  Machine  (SAM),  Chaotic  system,  Nonlinear  dynamics, 
Integrability,  Experimental  SAM  apparatus. 
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1.  Introduction 


This  paper  deals  with  the  Swinging  Atwood  Machine  (SAM),  a  non-linear  two- 
degrees-of-freedom  system  derived  from  the  well-known  simple  Atwood  machine.  The 
latter  was  devised  in  1784  by  George  Atwood,  a  London  Physics  lecturer  who  built  his 
own  apparata  as  a  means  of  practical  illustration,  in  order  to  experimentally  demonstrate 
the  uniformly  accelerated  motion  of  a  system  falling  under  the  earth  gravity  held  g  with 
mass  dependence  [5j.  In  Atwood’s  original  machine,  two  masses  are  mechanically  linked 
by  an  inextensible  thread  wound  round  a  pulley.  In  SAM ,  one  of  the  masses  (to)  is 
allowed  to  swing  in  a  plane  while  the  other  mass  (M)  plays  the  role  of  a  counterweight; 
hence  SAM  can  be  seen  as  a  parametric  pendulum  whose  length  is  varying  as  a  function 
of  the  parameter  g  =  M/m. 


For  about  twenty-five  years,  many  studies  have  been  conducted  concerning  the  me¬ 
chanical  behaviour  of  SAM.  Said  studies  were  conducted  exclusively  on  a  simplified  model 
for  SAM  neglecting  any  influence  from  a  massive  set  of  pulleys.  Through  numerical  inves¬ 
tigations,  [301]  inferred  the  pulley-less  SAM  to  be  an  extremely  intricate  system  exhibiting 
significant  changes  in  the  qualitative  behaviour  of  trajectories,  depending  on  /.t-values. 
Assuming  g  >  1,  motion  is  limited  in  space  and  two  types  of  trajectories  can  be  dis¬ 
tinguished  owing  to  the  initial  conditions:  singular  ones  for  which  pendulum  length  is 
initially  zero,  and  non-singular  ones  where  the  pendulum  is  initially  released  from  rest 
with  a  non-zero  length.  For  the  former,  it  appears  that  g  =  3  is  a  particular  condition 


>endulum  length  becomes 
29j.  The  latter  is  divided 


corresponding  to  terminating  trajectories,  i.e.  those  for  which 
zero  after  a  given  duration,  regardless  of  the  initial  conditions 
into  periodic,  quasi-periodic  and  what  can  be  conjectured  to  be  ergodic  trajectories  in 
some  domain.  SAM  without  massive  pulleys  was  also  studied  by  means  of  Poincare  sec¬ 
tions  wherein  chaotic  dynamic  behaviour  becomes  prominent  as  g  is  increased  [27|.  An 
interesting  and  surprising  result  is  the  integrability  of  the  pulley-less  SAM  for  g  =  3,  a 
conclusion  which  is  also  supported  using  Hamilton- Jacobi  theory  0]  and  Noether  sym¬ 
metries  0.  For  g  >  3,  [9]]  proved  that  SAM  without  massive  pulleys  is  not  integrable, 
contrary  to  what  was  speculated  by  0.  The  belonging  of  g  =  M/m  to  a  special  set  of 
parameters  {gp  :  p  G  Z}  was  established  as  a  necessary  condition  for  integrability;  this 
result  was  proven  independently  in  [§],  35 ]  and  Q,  and  is  proven  in  Remark  17. 11 21  of 


the  present  paper  as  well.  Moreover,  unbounded  trajectories  (g  <  1)  have  been  studied 
via  energetic  considerations  [31,];  0  identified  and  classified  all  periodic  trajectories  in 
the  pulley-less  SAM  for  g  =  3.  Finally,  a  very  recent  result  co-written  by  one  of  the 
authors  of  the  present  article  (0)  proved  the  non-integrability  of  this  pulley-less  model 
for  SAM  for  the  exceptional  values  {gp  :  p  G  Z };  this  had  been  an  open  problem,  at  least, 
since  Q  explicitly  tackled  the  issue  for  the  first  time.  It  is  worth  noting  that  all  of  these 
studies  are  theoretical,  albeit  for  the  most  part  strongly  supported  by  massive  numerical 
simulations. 


In  this  paper,  we  intend  to  describe  a  useful  physical  construction  of  SAM  in  detail, 
as  well  as  present  further  experimental  and  theoretical  results.  In  addition,  a  theoretical 
premise  is  introduced  which  stands  as  a  novelty  all  its  own:  to  wit,  as  suggested  by 
experiments,  pulleys  are  no  longer  neglected,  in  order  to  take  account  of  non-zero  radii 
and  rotation  around  their  axes  of  revolution.  When  dealing  with  7V-degree-of-freedom 
non-linear  systems,  the  modern  researcher’s  tendency  to  restrict  adjectives  such  as  “com- 
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plex”  to  N  3>  1  should  not  divert  us  from  the  fact  that  even  the  dynamics  of  apparently 
elementary  cases  such  as  N  =  2  are  often  very  difficult  to  determine  (0),  and  thus  a 
source  of  interest  in  their  own  right.  As  shown  in  this  paper,  such  is  the  case  for  SAM. 
A  schematic  representation  of  SAM  is  featured  in  Section  [5]  partly  aimed  at  the  deriva¬ 
tion  of  the  equations  of  motion  with  pulleys  in  Section  [3]  The  constructed  apparatus 
is  then  described  in  detail  in  Section  01  and  some  experimental  results  are  presented  in 
Section  [5j  A  comparison  with  the  theoretical  model  is  performed  in  Section  [6]  through 
numerical  simulations  of  the  general  equation  of  motion  obtained  in  Section  0J  Section  [3 
is  devoted  to  prove  the  non-integrability  for  SAM.  The  rigorous  proof  shown  therein  is 
requisite  to  establish  definitively  the  non-integrability  of  SAM,  hence  incumbent  upon 
any  proper  completion  of  experimental  and  numerical  results.  Indeed,  although  for  some 
systems  non-integrability  is  somehow  suggested  by  a  thorough  Poincare  section  anal¬ 
ysis,  chaotic  zone  detection  may  require  extremely  careful  numerics  and  will  at  times 
become  laborious  -  we  might  call  this  shy  chaos.  Furthermore  the  lack  of  integrability 
of  some  systems  cannot  be  discovered  by  looking  at  the  real  phase  space.  There  exist 
non-integrable  systems  without  any  recurrent  motion  in  the  real  phase  space  and  such 
that  chaos  is  confined  away  from  the  real  domain.  For  an  example,  see  [IS].  Finally, 
Section  [8]  concludes  with  perspectives  on  further  experiments  and  some  comments  on 
results  concerning  integrability  of  the  pulley-less  case. 

[30j  suggested  a  SAM  physical  demonstration  model  using  a  vertically  mounted  air 
table,  and  alleged  a  successful  experimental  demonstration  of  the  system’s  motions.  How¬ 
ever,  there  is  no  experimental  result  in  the  aforementioned  reference  and  its  proposed 
model  for  SAM  is  unequivocally  not  as  close  to  the  theoretical  system  as  the  model 
described  herein  (cf.  Section^-  To  our  present  knowledge,  detailed  experimental  studies 
of  SAM ,  let  alone  comparisons  of  any  such  experiments  with  the  theory,  do  not  exist  in 
the  literature  prior  to  our  work. 

Therefore,  our  work  arguably  completes  the  above  theoretical  and  experimental  re¬ 
search  on  SAM ,  and,  at  the  same  time,  opens  new  problems  and  sets  a  starting  point  for 
further  experimental  and  theoretical  studies. 

2.  Schematic  representation  of  SAM 

SAM  is  represented  by  the  system  S  sketched  in  Figure  [T]  and  consisting  of: 

•  a  pendulum,  considered  as  a  material  point  A  of  mass  to, 

•  a  counterweight,  considered  as  a  material  point  B  of  mass  M, 

•  a  thread  of  length  L  linking  A  and  B , 

•  two  pulleys  V  and  V'  of  radius  R,  distant  from  one  another  by  a  predetermined 
distance  D. 

S  is  studied  relative  to  the  Galilean  laboratory  frame  1Z  =  (O  ex  ey  ez)  whose  origin 
O  is  chosen  to  correspond  to  the  centre  of  the  pendulum  pulley.  Axis  Ox  corresponds 
to  the  pulley  revolution  axis;  axis  Oy  is  the  horizontal  direction  defined  by  the  pulley 
centres  (O  and  O'),  and  oriented  from  O  toward  O';  finally,  Oz  is  chosen  to  correspond, 
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Q 


B  (M) 


Figure  1:  Schematisation  of  the  system  S  =  {pendulum-thread-pulleys-counterweight}  representing 
SAM.  The  angle  cp  locates  the  material  point  N'  of  the  pulley  V . 


for  the  sake  of  convenience,  to  the  downward  direction  of  the  local  earth  gravity  field  g 
(vertical). 

Pendulum  A  is  characterized  by  its  variable  length  r  =  QA,  Q  being  the  geometrical 
point  where  the  thread  departs  from  the  pulley,  and  by  the  angle  9  formed  by  QA  and 
the  downward  vertical.  Note  that  9  as  represented  in  Figure  [T]  is  a  positive  angle. 

Vertical  motion  of  the  counterweight  B  is  described  by  its  coordinate  zb,  which  can 
be  related  to  the  angular  position  ip  of  any  point  N  on  the  pulleys,  provided  the  thread 
does  not  slip  on  the  pulley  (in  Figure  [1]  for  the  sake  of  clarity,  N  is  drawn  on  the  pulley 
V'  associated  to  B  and  thus  labelled  N').  Indeed,  under  this  assumption,  when  B  is 
falling  down,  pulleys  are  able  to  rotate  in  such  a  way  that  the  velocity  of  any  point  of 
the  pulleys  (for  instance  N' )  is  equal  to  the  velocity  of  B.  Hence: 

zB  =  Rp  (1) 

Note  the  difference  between  the  rotation  angle  p  of  the  pulleys  and  9:  the  former  defines 
the  location  of  any  material  point  on  a  pulley,  whereas  the  latter  defines  the  angular 
position  of  A,  as  well  as  that  of  the  geometrical  point  of  contact  Q.  This  subtlety  is  due 
to  the  necessary  mechanical  description  of  contact  in  terms  of  three  points  [2ll]: 

•  the  geometrical  point  of  contact  Q, 

•  the  point  N  of  the  pulley  V  and  corresponding  to  Q  at  time  t, 

•  the  point  K  of  the  thread  corresponding  to  Q  at  time  t. 

Figure [2] displays  the  configuration.  A  similar  problem  can  be  found  in  [22I]. 
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Figure  2:  Description  of  the  contact  between  a  pulley  and  the  thread  with  three  points:  K ,  Q,  and  N . 
At  any  time  t,  both  K  and  N  touch  each  other  at  Q.  In  the  figure,  these  points  are  separated  for  the 
sake  of  clarity. 


A  physical  way  of  understanding  the  difference  between  6  and  <p  is  to  imagine  the 
following  situation.  At  initial  time,  assume  that  Q  and  N  are  superposed:  0  =  Oq  and 
<p  =  tpo  (Figure Eli)-  If  0  is  fixed  and  B  heads  downwards  with  velocity  Zb,  the  absence 
of  slippage  of  the  thread  on  the  pulleys  implies  that  they  rotate  with  angular  velocity  p 
given  by  ©,  meaning  N  is  moving  while  Q  remains  fixed,  and  r  evolves  from  r o  to  r;  at 
final  time,  ip  ^  po  (Figure  EJd). 


Figure  3:  Evolution  of  point  N  of  the  pulley  V  when  6  is  kept  constant.  The  geometrical  point  of  contact 
Q  is  immobile  and  only  r  evolves. 


Since  L  is  supposed  to  be  constant,  Zb  is  directly  related  to  r,  D,  and  the  lengths 
QP  and  P'Q'  corresponding  to  regions  where  the  thread  and  the  pulleys  keep  contact. 


Because 


QP  =  R  (f  ~  °)  ’ 
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one  has,  precisely, 

L  =  r  +  D  +  7 tR  —  R9  +  zB  ■ 
Since  0  =  r  —  R9  +  zB,  it  follows  that 

z'b  RO  —  t 


R 


(2) 

(3) 


Hence,  S  is  a  system  with  two  degrees  of  freedom,  for  instance  9  and  r. 

3.  Equations  of  motion  of  SAM  with  pulleys 

3.1.  Equations  of  motion 

Let  us  determine  the  equations  of  motion  for  SAM  by  taking  into  account  the  pulleys, 
as  opposed  to  what  has  been  assumed  in  previous  theoretical  studies 


icount  tne  pulleys. 

liiii 


I30L 1311. 1351].  Indeed,  their  non-zero  radii  imply  a  likely  change  in  position  for  Q ,  and  the 
equally  likely  rotation  of  V  and  V'  around  their  respective  revolution  axes  apparently 
deems  their  inertial  momentum  Ip  a  significant  dynamic  parameter.  Observations  will 
confirm  this  -  see  Section  [5] 

Lagrange’s  formalism  is  used  to  derive  these  equations.  The  kinetic  energy  of  the 
system  is  expressed  by: 

£k  =  ^ mv\  +  ^ Mv2b  +  2  Qw 

The  first  term  sums  up  the  contribution  by  pendulum  A,  the  second  is  relative  to  the 
counterweight  B  and  the  third  one  corresponds  to  the  rotation  of  the  two  pulleys.  We 
have: 

AO  A  AOQ  AQA 

VA-^r~^r  +  ^T’ 

where  OQ  =  —Rcos9ey  —  Rsin9ez  and  QA  =  —rsin9ey  +  rcos9ez,  hence  in  the 
Cartesian  base  (O  ex  ey  ez)  we  can  write 


/ 


vA  = 


0 


\ 


R9  sin  9  —  r9  cos  9  —  r  sin  9 
y  —  R9  cos  9  —  rOsin.9  +  f  cos  9  ) 
Similarly,  vB  =  AOB/At  with  OB  =  zB  ez.  Using  ©,  one  gets 

vB  =  zBez  =  ( R9  -  f)  ez. 

Finally,  using  m  and  introducing  the  effective  total  mass  of  the  system 

2  Iv 

Mt=M  +  m  +  -jp, 


we  get 

This  expression  is  similar  to  that  obtained  when  neglecting  the  pulleys,  except  that: 


£k  =  l Mt{R6-f)2  +  ^mr26 2. 
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•  the  total  mass  is  now  different  from  M  +  to  by  the  term  2Ip/R2  conveying  the 
rotation  of  the  pulleys; 

•  the  counterweight  influences  pendulum  A  through  its  length  r  and  the  winding  of 
the  rope  on  the  associated  pulley.  The  latter  influence  is  considered  in  the  term 

Re. 

Potential  energy  is  only  due  to  the  Earth’s  local  gravity  field.  Dropping  an  irrelevant 
additional  constant  term,  we  have: 

£Pi9  =  —mg  ■  OA  —  Mg  ■  OB  ■  —mgzA  —  Mgzs, 


so  that 


£p>g  =  mg(R.  sin  9  —  r  cos  9)  +  Mg(r  —  R9). 

The  Lagrangian  L(r,  9,  r,  9)  =  £k  —  £p,g  of  the  system  is  thus: 

L(r,  9,  r,  9)  =  —Mt(R9  —  f)2  +  —mr292  —  gr(M  —  m  cos  9)  —  gR(m  sin  9  —  M9 ), 

from  which  one  deduces  the  conjugate  momenta  pr  and  pg  associated  to  r  and  9  respec¬ 
tively: 


pr  =  ~qT  =  —Mt(R9  -  r) 

(9t 

pg  =  — r  =  MtR(R9  —  r)  +  mr29  =  —Rpr  +  mr29 

89 

The  SAM  Hamiltonian  is,  in  this  case,  by  definition: 


'H  —  £k  +  £p,g  1 

or  else,  expressed  in  terms  of  pr  and  pg : 


H(r,9,pr,pg)  = 


1 


p2  (pg  +  Rpr)2 


Mt 


+  gr(M  —  to  cos  9)  +  gR(m  sin  9  —  AI9).  (4) 


Equations  of  motion  follow  from  the  Hamilton’s  equations: 


8H  ,  .  OH 

P’  =  -  37  »  =  -gp 

yielding 

(  pt(r  —  R9)  =  r92  +  g(cos9  —  p) 

\  r9  =  — 2r9  +  R92  —  g  sin  9 

with  /i  =  M/to  and  gt  =  Mt/m  =  1  +  g  +  (2 Ip/mR2). 


(5) 


9 


3.2.  A  more  physical  way  to  obtain  equations  of  motion 

An  alternative  method  can  be  used  in  order  to  derive  equations  of  motion  ([5]).  The 
second  of  these  is  obtained  by  applying  the  angular  momentum  theorem  at  the  mobile 
point  Q  in  order  to  cancel  reaction  force  at  this  contact  point  21] : 


AL 


Q 


At 


+  vq  x  m,VA  =  QA  x  mg , 


where  Lq  =  QA  x  mvA  =  mr29  ex  and 


RO 

'«  =  aqaq- 


v*  =  ikOQ+iAQA 


We  obtain  rO  =  —2 f6  +  R9 2  —  gsvnO.  Taking  this  into  account,  the  first  equation  in  © 
comes  from  the  conservation  of  mechanical  energy  in  SAM : 

£m  =  ^ Mt(R9  —  r )2  +  ^ mr29 2  +  mg(Rs'm9  —  rcos9)  —  Mg(R9  —  r)  =  C 

C  being  a  real  constant,  after  derivation  with  respect  to  t. 


3.3.  Comparison  with  previous  studies 

Without  any  pulley  influence,  i.e.  pulley  inertial  momentum  Ip  =  0  and  pulley  radius 
R  =  0,  we  recover  the  equations  obtained  by  [8(ij]: 


f  (1  +  A*)f_  =  r92  +  g{cos9  -  g)  ,  » 

\  r9  =  —2r9  —  gsm9 

Obviously,  if  there  is  no  oscillatory  motion  (9  =  0),  the  well-known  simple  Atwood 
machine  [5[  is  recovered: 

Mtr  =  g{m  —  M )  (7) 


4.  Description  of  the  experimental  apparatus 

A  physical  prototype  for  SAM  has  been  built  using  two  identical  pulleys,  a  nylon 
thread,  a  brass  ball  as  a  pendulum  and  a  set  of  different  hook  masses  acting  as  counter¬ 
weights.  The  pendulum  and  the  chosen  counterweight  are  linked  together  by  the  nylon 
thread  placed  around  the  pulleys.  A  photo  of  SAM  is  displayed  in  Figured 

4.1.  About  the  pendulum  and  the  counterweight 

Each  mass  in  the  experimental  device  has  been  measured  with  a  precision  scale  of 
0.01  g  of  accuracy.  The  pendulum  is  a  brass  ball  with  a  30  mm  diameter  and  a  mass 
m  =  118.36  g.  The  picture  of  the  pendulum  in  Figure  [5]i  also  exhibits  a  paper  clip  and 
the  nylon  thread,  the  latter  being  solidly  tied  to  the  brass  ball  and  the  paper  clip. 

The  paper  clip  is  secured  to  the  hook  of  the  chosen  counterweight,  in  turn  picked  out 
from  nine  hook  masses  whose  measured  values  are  M  =  10.01  g,  20.02 g(x2),  50.05 g, 
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Figure  4:  Photo  of  the  Swinging  Atwood  Machine  (SAM):  a  pendulum  (on  the  left)  and  a  counterweight 
(on  the  right),  linked  together  by  a  nylon  thread. 


100.10 g,  200.22 g(x2),  500.51  g,  and  1000. 10 g  (Figure [5Jd).  The  relative  difference  be¬ 
tween  these  values  and  those  engraved  in  each  hook  mass  is  0.1%;  thus,  with  respect 
to  the  orders  of  magnitude  of  the  different  masses  involved  in  the  experimental  device, 
this  difference  can  be  neglected.  The  values  considered  are  therefore  presumed  to  be 
those  indicated  on  the  hook  mass  themselves,  namely  M  =  10 g,  20g(x2),  50 g,  100 g, 
200g(x2),  500 g,  and  1000 g.  Henceforth,  and  for  the  sake  of  linguistic  simplicity,  these 
hook  masses  will  be  called  “counterweights” ,  although  weight  and  mass  are  different  no¬ 
tions,  however  related.  This  set  enables  varying  the  counterweight  mass  from  10  g  (one 
mass)  to  2  100  g  (addition  of  all  the  masses)  with  a  step  of  10  g  by  hooking  several  masses 
together.  Among  these  hook  masses,  one  is  hung  on  the  nylon  thread  by  means  of  the 
paper  clip,  whose  measured  mass  is  0.37  g.  It  is  interesting  to  note  that,  by  a  fortunate 
coincidence,  the  mass  of  the  paper  clip  is  equal,  with  a  0.01  g  difference,  to  0.36  g,  i.e. 
the  mass  of  the  brass  ball  minus  118  g.  Therefore,  the  mass  of  the  brass  ball  can  be  taken 
as  equal  to  118  g  and  the  mass  of  the  paper  clip  can  be  ignored.  Finally,  we  get,  for  the 
pendulum  and  the  counterweight,  respectively:  m  =  118  g  and  10  g  <  M  <  2100  g. 

4-2.  About  the  nylon  thread  and  the  pulleys 

The  thread  (Figure  [5ji)  ensures  a  mechanical  coupling  between  the  pendulum  and 
the  counterweight  through  the  two  pulleys.  The  length  of  the  thread  is  about  one  meter 
and  its  measured  mass  of  0.10  g  is  negligible  compared  to  the  other  masses  involved.  In 
addition,  the  nylon  thread  is  assumed  inextensible.  During  experimentation,  no  thread 
breaking  has  been  reported. 

Pulleys  used  are  shown  in  photos  of  Figure  [Gj  They  are  made  up  of  two  parts: 
an  internal,  immobile  one  bound  to  the  revolution  axis,  and  a  mobile,  external  one 
liable  to  rotate  around  this  axis.  These  two  pulley  components  are  uncoupled  through  a 
ball  bearing  which,  moreover,  reduces  mechanical  energy  dissipation  by  friction.  Pulley 
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Figure  5:  (a)  Brass  ball  (pendulum),  nylon  thread  and  paper  clip  (b)  Set  of  hook  masses  used  in  the 
experiment  (counterweight).  From  left  to  right:  10 g,  20 g,  20 g,  50 g,  100  g,  200 g,  200  g,  500 g,  1000 g. 


radius  is  R  =  2.5  cm  and  that  of  the  motionless  part  is  1  cm.  Pulley  V,  associated  to  the 
pendulum,  has  been  modified  in  order  to  make  its  groove  deeper.  Indeed,  during  the  first 
experimentations  we  observed  that  the  thread  could  rapidly  exit  the  groove  because  of 
the  pendulum  motion.  To  avoid  this,  which  could  by  the  way  be  dangerous,  two  metallic 
plates  were  added  and  fixed  to  the  pulley  in  order  to  increase  by  1  cm  the  depth  of  the 
groove  (Figures  [6ji  and  [6jo) .  It  is  worth  noting  that  the  plates  are  fixed  to  the  immobile 
part  of  the  pulley  and  are  in  no  way  in  contact  with  the  mobile  one.  The  motion  of  the 
latter  one  is  thus  not  affected  by  such  a  modification:  hence,  from  a  mechanical  point  of 
view,  the  resulting  pulley  is  identical  to  the  original  one. 

f.3.  Strengthening  of  the  apparatus 

Figure  2]  also  features  two  horizontal  metallic  rods  binding  the  two  feet  of  SAM. 
Their  role  is  to  reinforce  the  machine.  Indeed,  due  to  considerable  stress  involved  in  the 
pendulum  and  counterweight  motions,  the  orientation  of  the  two  pulleys,  as  well  as  the 
arbitrary  distance  D  =  57.5  cm  between  them,  can  change;  thus,  it  could  be  dangerous 
not  to  strengthen  the  whole  device.  The  lower  metallic  rod  is  solidly  fixed  to  the  feet 
while  the  upper  one  is  solidly  fixed  to  the  revolution  axes  of  the  pulleys.  Consequently, 
we  ensured  a  constant  distance  between  the  pulleys  whose  axes  keep  a  constant  direction; 
in  virtue  of  such  a  construction,  SAM  is  solid  and  operational. 


5.  Motion  of  the  pendulum:  experimental  results 


5.1.  Experimental  measure  of  Ip 

The  presence  of  the  inertial  momentum  Ip  of  a  pulley  in  (J5])  renders  an  experimental 
determination  thereof  necessary.  This  measure  was  made  using  a  simple  Atwood  machine 
where  the  heavier  mass  (M  =  130  g)  fell  down  from  a  convenient  height  h  =  1.66  m;  the 
lighter  mass  (the  brass  ball)  being  m  =  118  g.  Using  equation  J7|),  one  gets: 


IP  =  R2 


( M  —  m)g(/S.ty 


(M  +  in) 


2  h 
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Figure  6:  Photos  of  the  two  pulleys,  (a)  Pulley  V ,  front  view:  a  metallic  plate  is  added  which  hides 
the  internal  part  and  the  ball  bearing  of  the  original  pulley,  (b)  Pulley  P,  top  view:  the  groove  is  made 
deeper  by  1  cm  with  the  additional  plate;  the  original  pulley  is  easily  recognizable  between  the  two  plates, 
(c)  Pulley  V' ,  front  view:  the  ball  bearing  uncouples  the  immobile  internal  part  and  the  mobile  external 
one  (d)  Pulley  Px,  top  view:  photograph  of  the  axis  of  revolution  and  the  groove  of  the  pulley.  Moreover, 
the  two  pulleys  are  distant  enough  from  one  another  to  avoid  a  pendulum-counterweight  collision  during 
the  motion. 


where  At  is  the  fall  duration.  Through  a  set  of  ten  measures  with  a  chronometer  of 
0.01  s  of  accuracy,  the  mean  fall  duration  found  is  (At)  =  2.70  ±  0.01  s.  (Ip)  =  6.85  x 
10~6kg  •  nr2  ensues.  Concerning  the  value  of  Ip,  errors  are  due  to  the  measure  of  both 
At  and  the  positions  of  M  at  the  initial  and  final  times  -  that  is,  the  determination 
of  h.  Uncertainties  are  mainly  due  to  the  determination  of  the  final  time,  which  must 
correspond  to  the  falling  distance  h  as  precisely  as  possible;  initial  and  final  positions  are 
determined  with  an  error  of  0.1  cm,  which  compared  to  the  value  of  h  can  be  neglected. 
Consequently,  uncertainties  in  position  determinations  are  disregarded  and  the  error 
in  the  measure  of  Ip  can  be  reasonably  associated  to  the  uncertainty  in  At  (0.4  %). 
Thus,  precision  on  Ip  is  twice  that  of  At,  hence  about  1  %;  absolute  uncertainty  is  thus 
0.07  x  10  6  kg  •  m2.  Therefore,  we  can  write: 

Ip  =  (6.85  ±  0.07)  x  10"6  kg  •  m2  or  Ip  =  6.85  x  10~6  kg  •  m2  ±  1  %. 

5.2.  Experimental  results 

The  motion  of  the  pendulum  has  been  filmed  for  various  ^-values  and  initial  conditions 
(ro  ;  $o).  Then,  using  the  “Synchronie”  software  and  focusing  on  each  film  image  by 
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image,  a  pointer  enabled  us  to  pick  up  pendulum  positions  and  record  them.  Such  a 
process  is  necessarily  a  source  of  errors,  as  it  is  sometimes  difficult  to  locate  the  pendulum 
exactly,  especially  if  velocity  is  high.  The  errors  introduced  by  such  a  procedure  are  not 
simple  to  estimate.  However,  trajectories  have  been  correctly  recorded,  as  comparisons 
with  numerically-simulated  theoretical  results  will  show  (see  Section  [6]) . 

5.2.1.  Case  p.theo  =  3 

Since  previous  studies  focused  mainly  on  the  particular  and  theoretical  case  ptheo  =  3, 
this  was  the  first  one  we  experimentally  addressed.  In  fact,  the  masses  available  only 
allowed  us  to  approach  p,theo-  with  m  =  118  g  and  M  =  350  g  one  obtains  p.exp  = 
2.966,  which  is  the  closest  value  to  ^theo-  The  sampling  time  step  has  been  67  ms. 
Motion  has  been  researched  for  four  different  initial  conditions  (ro  ;  0o)=(O.649;  53.5), 
(0.710;  66.5),  (0.854;  68.3)  and  (0.867 ;  51.1);  ro  is  in  meters,  90  in  degrees.  The  motion 
of  the  pendulum  presents  the  same  pattern  and  characteristics  for  all  these  conditions, 
so  only  the  trajectory  for  the  first  initial  conditions  is  shown  in  Figure  [7Ji.  All  in  all, 
359  sampling  times  have  been  recorded.  The  pendulum  has  a  planar  revolving  trajectory 
around  the  pulley  and  presents  an  asymmetry  with  respect  to  the  vertical  direction. 
Note  that  the  pendulum  becomes  closer  and  closer  to  the  pulley  as  a  consequence  of 
dissipative  phenomena  and  is  bound  to  end  up  knocking  against  it.  Phenomena  qualifying 
as  dissipative  are,  to  our  knowledge:  the  friction  between  the  thread  and  the  pulley, 
the  air  friction  on  the  pendulum  and  the  counterweight  as  well  as  friction  inside  the 
ball  bearing.  Evolutions  of  the  length  of  the  pendulum  r  and  angle  9  are  displayed  in 
Figure  [7])  and  Figure  [T):  respectively.  The  asymmetry  of  the  trajectory  and  dissipation 
are  observable  in  the  evolution  of  r  since  this  variable  exhibits  different  minimal  and 
maximal  amplitudes  which  decrease  in  function  of  time  t. 

The  Fourier  analysis  (non-displayed)  of  the  data  shows  that,  for  the  behaviour  of 
9 ,  the  most  relevant  harmonic  is  the  constant  term,  followed  by  harmonics  7,  6  and 
8.  They  account,  respectively,  for  0.401,  0.311,  0.117,  and  0.091  of  the  total  variation. 
Of  course,  the  contributions  of  harmonics  6  and  8  are  due  to  leakage  of  the  7th  one. 
The  amplitude  of  the  6th  harmonic  is  larger  than  that  of  the  8th,  showing  that  the 
true  dominant  average  frequency  is  slightly  less  than  7  times  the  basic  frequency.  From 
the  number  of  data  and  the  sampling  step  time,  a  basic  frequency  of  0.042  Hz  follows. 
Hence,  an  average  value  of  the  dominant  frequency  can  be  estimated  to  be  equal  to 
0.294  Hz.  However,  as  is  clear  by  looking  at  the  maxima  of  the  angle,  the  frequency 
changes  with  time.  The  spacing  between  successive  maxima  takes  the  approximate  val¬ 
ues  3. 820s, 3. 909s, 3. 379s, 3. 311s, 3. 173s, 2. 780s,  and  2.643s.  Hence,  the  instantaneous 
frequency  changes  from  about  0.262  Hz  to  about  0.378  Hz.  The  explanation  is  simple: 
the  dissipation  reduces  the  energy  and  the  length  of  the  pendulum  becomes  shorter, 
increasing  the  frequency. 

For  the  radius,  the  major  contribution  comes  from  the  constant  term,  followed  by 
harmonics  13  and  7.  They  account  for  0.810,  0.037,  and  0.026  of  the  total  variation. 
Comparing  the  plot  of  r  as  a  function  of  time  with  that  of  9,  one  can  see  a  doubling 
in  the  number  of  maxima.  The  explanation  is  simple:  largest  maxima  occur  at  the  left 
part  of  the  plot  of  the  orbit.  Then,  a  minimum  is  reached  when  9  =  —  7r  (i.e. ,  upwards), 
followed  by  a  maximum  to  the  right,  and  a  new  minimum  at  9  =  —  7r  to  reach  a  larger 
maximum  on  the  left.  The  spacing  between  successive  larger  maxima  of  r  is  very  close  to 


14 


the  one  observed  for  9.  As  mentioned,  the  largest  non-constant  harmonics  are  the  13th 
and  the  7th,  not  in  a  2-to-l  ratio.  This  is  related  to  the  fact  that  the  “best”  estimate 
of  the  main  frequency  for  9  is  slightly  less  than  7  times  the  basic  frequency.  The  closest 
integer  to  the  double  would  be  13  rather  than  14. 

The  decrease  in  energy  can  be  calculated  as  follows.  From  the  experimental  data,  the 
values  of  r  and  9  can  be  computed.  For  that,  we  have  used  two  independent  methods. 
The  first  one  is  simply  numerical  differentiation  with  a  central  formula.  The  second 
one  aims  at  filtering  errors  in  the  data  as  well.  A  Discrete  Fourier  Transform  has  been 
computed  and  harmonics  up  to  order  32  have  been  retained.  Then,  it  is  possible  to  check 
that  the  reconstruction  agrees  quite  well  with  the  initial  data  and  one  can  compute  the 
values  of  r  and  9  using  these  Fourier  expansions.  To  prevent  leakage  due  to  the  fact  that 
the  data  at  the  ends  of  the  interval  are  quite  different  (Figs.  [7)3  and  c),  which  originates 
the  well-known  0(l/n)  decrease  of  the  order  of  magnitude  of  the  nth  harmonic,  different 
procedures  have  been  used,  but  the  results  are  essentially  the  same.  They  also  show  a 
reasonable  agreement  using  the  first  and  second  methods. 


Figure  7:  Case  p.eXp  =  2.966  with  initial  conditions  ro  =  0.649  m  and  Qo  =  53.5°.  (a)  Experimental 
positions  (black  crosses)  and  interpolated  pendulum  trajectory  (solid  line).  The  Cartesian  coordinates 
of  the  initial  position  are  y  =  —  0.54  m  and  z  =  0.37  m  and  those  of  the  final  position  are  y  =  0.02  m  and 
z  =  —0.05  m.  The  pulley  is  represented  by  the  circle  whose  centre  is  at  the  origin  of  coordinates  (0 ;  0). 
(b)  Experimental  positions  (black  crosses)  and  interpolated  curve  (solid  line)  for  the  evolution  of  r.  (c) 
Same  as  (b)  but  for  the  angle  6. 
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When  r  and  9  are  available,  one  can  compute  pr  and  pg  and  subsequently  the  value 
of  the  energy.  The  values  for  which  9  reaches  a  maximum  (i.e. ,  on  the  left  of  Figure^) 
are  shown  in  Figure  [3]  The  rate  of  decrease  of  the  energy  is  about  0.037  J  •  s-1. 
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Figure  8:  Values  of  the  energy  (crosses)  occurring  when  6  reaches  a  maximum  at  the  left  part  of  Figure[7}i. 
We  have  also  displayed  the  function  £rn  =  i.723  —  0.037 1. 


5.2.2.  Case  ptheo  =  1.5 

The  experimental  value  of  p  closest  to  ptheo  =  1.5  is  /ieXp  =  1-525,  given  by  M  =  180  g. 
Proceeding  as  in  the  above  case  allows  us  to  retrieve  the  experimental  trajectory  and 
the  evolution  of  the  degrees  of  freedom.  In  this  case,  time  step  is  40  ms.  Two  initial 
conditions  have  been  considered:  (ro ;  0o)=(O-484m;  87.0°)  and  (0.621m;  87.7°).  Since 
they  produce  the  same  dynamic  behaviour,  only  the  first  one  is  displayed  (Figures 
and[5t  respectively).  A  slight  asymmetric  trajectory  with  respect  to  the  vertical  direction 
and  an  evolution  of  9  close  to  periodic  with  a  period  around  1.1  —  1.2  s  can  be  observed. 
For  the  evolution  of  r,  asymmetry  and  slight  dissipation  are  also  observed. 

A  study  similar  to  that  of  ^texp  =  2.966  is  performed.  The  main  contribution  to  the 
Fourier  analysis  of  9  comes  from  harmonic  number  7  which  accounts  for  0.912  of  the 
total  variation  with  a  frequency  of  approximately  0.875  Hz,  although  a  better  value  for 
the  average  frequency  seems  to  be  0.89  Hz.  Harmonics  14  and  21  also  play  a  relevant 
role.  As  we  did  for  /reXp  =  2.966,  we  can  consider  the  spacing  between  successive  maxima 
of  9  which  takes  the  values  1.125  s,  1.119  s,  1.114  s,  1.110  s,  1.156  s,  1.133  s  and  1.063  s, 
showing  a  decreasing  trend  with  irregularities. 

For  r,  the  largest  harmonic  is  the  constant  term  which  accounts  for  0.971  of  the  signal. 
If  we  skip  this  term,  harmonics  7,  14  and  3  are  clearly  seen.  They  contribute  to  0.30, 
0.29,  and  0.14  of  the  signal  minus  the  constant  part. 
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(a) 


(b) 


t(s) 

(c) 


Figure  9:  Same  as  Figure[7]but  for  ^eXp  =  1.525  and  initial  conditions  ro  =  0.484m  and  0q  =  87.0°.  In 
(a),  the  Cartesian  coordinates  of  the  initial  position  are  y  =  —0.485  m  and  z  =  0  ;  for  the  final  position, 
one  has  y  =  —0.330  m  and  z  =  0.195  m.  The  origin  (0;  0)  corresponds  to  the  centre  of  the  pulley  (not 
represented) . 
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The  decrease  of  the  energy  as  a  function  of  time  has  been  displayed  in  Figure  fTT)l  this 
time  using  the  values  of  the  energy  computed  at  the  minima  of  9,  on  the  right  part  of 
Figure  [HJi-  Now  the  rate  of  decrease  is  about  0.024  J  s-1.  In  this  case,  using  the  filtered 
Fourier  methods  gives  better  results,  because  of  the  large  changes  in  position  with  a  time 
step  of  40  ms. 


t(s) 


Figure  10:  Values  of  the  energy  (crosses)  occurring  when  6  reaches  a  minimum  at  the  right  part  of 
Figure  j9ji.  We  have  also  displayed  the  function  £m  =  0.751  —  0.024 1. 


5.2.3.  An  unbounded  case:  p.theo  =  1 

In  this  situation,  the  experimental  value  of  /. i  is  /reXp  =  1-017  and  three  initial  con¬ 
ditions  have  been  considered:  tq  =  0.120m,  0.263m  and  0.477m  for  9q  =  77.9°.  Again, 
trajectories  present  the  same  pattern,  so  only  one  is  shown  (Figure  lllhl.  They  are 
characterized  by  an  increase  in  r  and  9  oscillations  around  the  vertical  [9  =  0)  with  a 
decreasing  amplitude  of  9  (Figures  [TTJd  and  c  respectively).  In  Figure [TTb.  r  appears  to 
approach  a  linear  increase  in  time:  r  =  vzt  where  vz  is  the  velocity  along  the  vertical. 
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(a) 


y  (m) 
(b) 


(c) 


Figure  11:  Same  as  Figure[7]but  for  y,ex p  =  1-017  and  initial  conditions  ro  =  0.120m  and  6q  =  77.9°.  In 
(a),  the  Cartesian  coordinates  of  the  initial  position  are  y  =  —0.122  m  and  z  =  0  ;  for  the  final  position, 
one  has  y  =  —0.116  m  and  z  =  1.167  m.  The  origin  (0;  0)  corresponds  to  the  centre  of  the  pulley  (not 
represented) . 
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6.  Numerical  solution  of  SAM  equations  of  motion 

6.1.  Theoretical  trajectories 

Equations  of  motion  ([5])  have  been  numerically  integrated  for  the  same  initial  condi¬ 
tions  and  values  of  the  parameter  /i  as  above  in  order  to  compare  the  theoretical  trajecto¬ 
ries,  displayed  in  Figure  fl^l  with  the  experimental  ones.  Computation  of  the  theoretical 
trajectories  has  been  performed  by  using  both  a  Taylor  integration  method  and  a  vari¬ 
ety  of  Runge-Kutta  methods  of  different  orders  with  step-size  control  to  integrate  the 
equations  of  motion. 

From  a  general  point  of  view,  the  theoretical  trajectories  seem  quite  similar  to  the 
experimental  ones.  However,  some  slight  differences  can  be  detected.  Firstly,  it  is  obvious 
that  dissipative  phenomena,  though  experimentally  reduced,  play  a  non-negligible  role 
since  convergence  of  the  pendulum  towards  the  pulley  for  the  first  case  (Fig.  [7b),  decrease 
of  r  for  the  second  one  (Fig.  (5Jd)  and  relatively  slow  increasing  of  r  for  the  third  one 
(Fig.  mb')  are  clearly  associated  to  energy  dissipation.  Friction  will  be  studied  a  bit 
further  in  Subsection 03  and  much  more  in  future  works.  Secondly,  it  must  be  said  that 
for  the  case  /rexp  =  2.966  the  pendulum  touches  the  nylon  thread  at  each  revolution,  an 
effect  which  could  be  included  into  the  equations  of  motion  through  a  dissipative  term; 
however,  it  seems  quite  difficult  to  introduce  such  an  effect  in  a  realistic  manner. 

An  important  point  concerns  the  comparison  between  these  trajectories  and  those  of 
i(3ij].  For  the  first  case  for  instance,  Tufillaro’s  trajectories  are  symmetrical  with  respect 
to  the  vertical  direction,  as  opposed  to  the  above  ones  (Figures  [7]i  and  fl2k').  Clearly, 
this  asymmetry  is  due  to  the  influence  of  the  pulleys. 

6.2.  Influence  of  the  pulleys  on  the  motion 

Pulleys  can  influence  the  motion  of  SAM  through  their  dimension  (since  radius  R  0) 
and  their  rotation  (since  Ip  0). 

6.2.1.  Influence  of  the  radius:  pexp  =  2.966 

Figure [13]  sketches  the  trajectory  for  /rexp  =  2.966  when  Ip  =  0  for  different  values  of 
increasing  i?;  i.e.  the  pendulum  pulley  V  has  a  non-negligible  radius  but  pulleys  are  not 
allowed  to  rotate.  The  first  figure  corresponds  to  the  symmetrical  Tufillaro  trajectory 
( R  =  0,  Ip  =  0).  Ostensibly,  the  larger  R  is,  the  more  significant  the  asymmetry  becomes; 
for  R  =  5  cm,  the  pendulum  hits  the  pulley  before  completing  one  revolution. 

6.2.2.  Influence  of  the  inertial  momentum:  fj.exp  =  2.966 

Figure  [14]  sketches  the  trajectory  for  /rexp  =  2.966  for  increasing  values  of  Ip  with 
a  value  of  pulleys  radius  taken  to  the  real  value  R  =  2.5  cm.  The  trajectory  of  the 
pendulum  is  visibly  modified:  it  describes  more  irregular  trajectories  and  fills  more  space 
as  Ip  increases. 


z  (m)  2  (m)  z  (m) 
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(a) 


(b) 


(c) 


Figure  12:  Numerical  trajectories  obtained  by  solving  equations  of  motion  ©.  (a)  //exp  =  2.966  and 
for  the  same  initial  conditions  as  Figure  [7^.  (b)  //eXp  =  1.525  and  for  the  same  initial  conditions  as 
Figure  [§Ji.  (c)  //eXp  =  1.017  and  for  the  same  initial  conditions  as  Figure  ITOi. 
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(a) 


(b) 


Figure  13:  Theoretical  trajectories  of  the  pendulum  with  fixed  pulleys  of  different  radii  R  for  /xexp  =  2.966 
and  initial  condition  (ro  ;  0q)= (0.649m;  53.5°).  (a)  R  =  0,  (b)  R  =  2.5cm,  (c)  R  =  4cm  and  (d) 
R  =  5  cm.  In  the  latter  subplot,  the  pulley  has  been  represented  by  the  circle  whose  centre  is  at  the 
origin  of  coordinates  in  order  to  see  the  collision  between  the  pendulum  and  the  pulley. 

(a)  (b) 


Figure  14:  Theoretical  trajectories  of  the  pendulum  with  pulleys  of  radius  R  =  2.5  cm  and  different 
inertial  momenta  Ip  for  /xeXp  =  2.966  and  initial  condition  (ro  ;  #o)=(0.649m;  55.7°).  (a)  Ip  =  0;  this 
subplot  which  is  the  same  as  Figure  Il3b  has  been  repeated  for  the  sake  of  clarity,  (b)  Ip  =  6.85  x 
10-6  kg  •  m2,  (c)  Ip  =  13.70  x  10-6  kg  •  m2  and  (d)  Ip  =  27.40  x  10-6  kg  •  m2. 
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(a) 


(c) 


(b) 


(d) 


Figure  15:  Theoretical  trajectories  of  the  pendulum  with  fixed  pulleys  of  different  radii  R  for  /xex p  =  1-525 


and  initial  condition  (ro  ;  0q)— (0.484  m;  87.0°). 
R  =  10  cm. 


(a) 


(c) 


(a)  R  =  0,  (b)  R  =  2.5  cm,  (c)  R  =  5  cm  and  (d) 


(b) 


(d) 


Figure  16:  Theoretical  trajectories  of  the  pendulum  with  pulleys  of  radius  R  =  2.5  cm  and  different 
inertial  momenta  Ip  for  /xeXp  =  1.525  and  initial  condition  (ro  ;  #o)=(0.484m;  87.0°).  (a)  Ip  =  0;  this 
subplot  which  is  the  same  as  Figure  Il5b  has  been  repeated  for  the  sake  of  clarity,  (b)  Ip  =  6.85  x 
10-6  kg  •  m2,  (c)  Ip  =  13.70  x  10-6  kg  •  m2  and  (d)  Ip  =  54.80  x  10-6  kg  •  m2. 
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6.2.3.  Influence  of  the  pulleys:  pex p  =  1.525 

In  this  case,  results  of  variations  in  R  with  Ip  =  0  (Figure  fl5l)  and  Ip  variations  for 
R  =  2.5  cm  ('Figure|  161)  are  similar  in  that  trajectories  are  modified  more  visibly  as  R  and 
Ip  increase.  However,  the  influence  of  Ip  and  R  on  SAM  motion  depends  on  the  value 
of  n  considered.  If  R  increases,  the  pulleys  being  fixed,  the  brass  ball  ends  up  hitting 
the  pulley  but  asymmetry  does  not  becomes  more  and  more  important.  If  Ip  increases, 
pulley  dimensions  being  fixed,  the  pendulum  evolves  in  a  much  more  limited  space. 

6.3.  Poincare  maps  and  rotation  number 

To  have  a  global  view  of  the  dynamics  of  SAM ,  we  have  computed  Poincare  maps 
Vm  on  suitable  Poincare  sections.  We  note  that  the  Hamiltonian  (HJ)  is  not  27r-periodic 
because  of  the  linear  term  in  9.  Be  that  as  it  may,  our  Poincare  section  E  defining  Vm 
is  given  by  the  coiling  of  6  through  multiples  of  27r,  with  9  >  0.  But,  contrarily  to  SAM 
without  pulleys  ([23]),  one  has  t°  distinguish  between  cuts  through  different  multiples  of 
27t.  Plus,  one  can  not  superimpose  the  different  sheets. 

Some  types  of  “escape”  can  occur.  The  main  source  thereof  is  r  going  to  zero.  Other 
relevant  sources  of  escape  are  r  increasing  too  much  or  \0\  becoming  too  large.  All  orbits 
leading  to  some  of  these  escapes  are  deleted. 

To  compare  with  the  experiments,  we  present  some  examples  for  M  =  350  g  and  M  = 
180  g.  As  levels  of  energy,  one  has  taken  the  values  corresponding  to  the  experiments  de¬ 
scribed  in [AlT] and [5221  that  is,  (tq  ;  0o)=(O.649m;  53.5°)  and  (ro  ;  0o)=(O-484m;  87.0°), 
respectively,  with  zero  initial  velocity.  Given  r,  pr  in  E  and  9  =  0,  the  value  of  pg  is 
recovered  from  the  energy  level.  Figure  [H]  shows  some  results.  All  these  massive  com¬ 
putations  use  Taylor  integration  methods,  in  order  to  ensure  a  very  good  conservation 
of  the  energy  (see,  e.g.,  |24j). 

The  top  left  plot  corresponds  to  M  =  350  g,  leaving  E  with  pr  >  0.  The  points  on 
E  with  pr  <  0  correspond  to  cuts  through  9  =  — 2-7T,  in  agreement  with  the  description 
of  motion  in  15.2.11  To  produce  the  Poincare  map,  we  first  computed  the  periodic  orbit 
as  a  fixed  point  in  ED  {pr  >  0}.  The  approximate  values  (r*  ;  p*)  of  the  fixed  point  are 
(0.332814  m ;  0.554330  kg  •  m  •  s-1).  This  periodic  orbit  can  also  be  obtained  by  starting 
the  motion  from  rest  at  (ro  ;  9q)  «  (0.61316  m;  64.032°),  not  too  far  from  the  values 
used  in  the  experiment.  The  curve  drawn  with  large  dots  around  the  periodic  orbit 
shows  the  iterates  of  Vm  corresponding  to  the  data  used  in  the  experiment.  It  is  clear 
that  the  theoretical,  non-dissipative,  motion  seems  to  be  in  a  2D  torus,  but  numerical 
computations  can  never  exclude  the  possibility  of  having  a  periodic  orbit  with  very  long 
period  or  a  tiny  chaotic  zone.  The  intersection  of  this  torus  with  9  =  — 27r  is  also  shown 
in  the  region  pr  <  0.  The  asymmetry  is  ostensible  in  the  plot,  and  is  due  to  the  effect  of 
the  pulleys.  To  produce  the  full  plot,  we  have  taken  initial  conditions  on  E  with  r  =  r* 
and  different  values  of  pr  starting  at  p*.  From  some  value  of  pr  onwards,  iterates  escape. 
We  do  not  exclude  the  presence  of  tiny  islands  outside  the  last  invariant  curve  shown. 

For  the  sake  of  completeness,  we  show  in  the  top  right  plot  the  Poincare  iterates 
leaving  E  in  the  region  {pr  <  0}.  Points  appearing  in  {pr  >  0}  are  on  the  sheet  {9  =  2n}. 
The  fixed  point  in  {pr  <  0}  is  approximately  (0.383367m;  -0.620020kg-  nr  •  s-1). 
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Figure  17:  Iterates  of  Poincare  maps  for  M  =  350  g  (top  row)  and  M  =  180  g  (bottom  plot)  on  the 
levels  of  energy  of  the  physical  experiments  described  in  the  text.  The  variables  (r ;  pr ),  as  usual  in  m 
and  kg  •  m  •  s— 1,  respectively,  are  shown  for  6  =  0  (mod  2n)  and  0  >  0.  On  the  top  left  plot,  the  points 
leave  the  section  with  pr  >  0  and  6  =  0,  and  those  with  pr  <  0  are  in  6  =  —2tv.  On  the  right  top  row, 
top  initial  points  are  in  6  =  0  with  pr  <  0,  and  those  showing  in  pr  >  0  have  6  =  2tt.  On  the  bottom 
plot,  all  intersections  occur  at  9  =  0.  The  thick  curves  are  the  intersections  of  the  theoretical  2D  tori 
corresponding  to  the  physical  experiments  with  the  section  E. 

The  bottom  plot  corresponds  to  M  =  180  g.  In  that  case,  only  intersections  having 
6  =  0  are  found.  The  fixed  point  is  (r*  ;  p*)  «  (0.356386  m;  0.008848kg-  m  •  s_1),  which 
can  also  be  obtained  leaving  from  (ro  ;  0O)  =  (0.53313  m;  78.596°),  again  not  too  far 
from  the  values  used  in  the  experiment.  As  before,  the  curve  drawn  with  large  points 
around  the  periodic  orbit  would  be  the  one  obtained  for  the  physical  experiment  without 
dissipation  and,  as  expected,  denotes  motion  on  a  2D  torus. 

A  useful  tool  to  understand  the  dynamics  of  Area-Preserving  Maps  and,  in  particular, 
Poincare  maps  such  as  the  ones  displayed,  is  the  rotation  number  p  for  the  map  restricted 
to  invariant  curves.  Despite  the  fact  that  the  rotation  number  still  exists  for  periodic 
orbits  of  Vm  and  for  the  eventual  islands  around  them  (thereupon  being  rational),  it 
is  not  defined,  in  general,  for  orbits  with  chaotic  dynamics.  The  method  used  for  the 
computation  is  topological  and  based  on  the  order  of  the  arguments  of  the  iterates  with 
respect  to  the  central  fixed  point  (r*  ;  p*)  of  Vm.  The  procedure  computes  two  estimates 
Pinf  and  psup  such  that  pmf  <  p  <  psup.  If  for  some  orbit  one  has  p-m{  >  psup  this  proves 
p  not  defined.  If  the  number  of  Pm  iterates  is  N,  the  typical  errors,  when  p  exists,  are 
O  (-V-2)  for  constant  type  rotation  numbers.  See  the  Appendix  in  [23J  for  details  and  a 
complete  analysis  of  the  error  depending  on  the  Diophantine  properties  of  p. 
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In  Figure  fl8l  we  show  results  corresponding  to  the  Poincare  maps  displayed  in  Figure 
E3  top  and  bottom.  The  computations  are  done  starting  at  initial  points  of  the  form 
(r*  ;  pr)  with  pr  =  p*  —  j A,  j  =  1,2,...  with  a  small  step  A.  As  successive  iterates  fall 
in  E  for  values  of  9  alternating  between  0  and  —  27r,  the  map  'P./,  has  been  used  instead 
of  Vm.  On  the  left  plot,  we  display  p  as  a  function  of  pr  for  M  =  350  g.  This  is  the  curve 
which  has  a  large  dot  near  the  upper  right  corner.  The  point  corresponds  to  pr  =  p*  and 
the  limit  rotation  number.  We  see  a  decreasing  rotation  number  when  pr  decreases  up  to 
a  value  pr  ss  0.143982  kg  •  m  •  s-1.  Beyond  that  point,  the  Vm  iterates  escape.  The  other 
curve,  also  shown  here  for  comparison,  corresponds  to  the  pulley-less  case  with  p  =  3 
and  will  be  described  later. 

In  fact,  what  seems  a  nice  curve  for  M  =  350  g  should  have,  generically,  a  “devil’s 
staircase”  structure.  That  is,  there  are  infinitely  many  intervals  in  which  p  £  Q  and  it  is 
constant.  They  correspond  to  islands  around  elliptic  fixed  points.  Some  of  these  rational 
values,  such  as  1/6,  1/7,  1/9,  3/20,  ...  (or  resonances)  have  been  detected.  But  they  are 
very  narrow.  As  an  example,  the  inset  in  the  left  plot  of  Figure  [TS]  shows  the  behaviour 
of  p  in  an  interval  of  pr  whose  width  is  lower  than  10_3kg  •  m  •  s^1  and  p  around  1/9, 
which  illustrates  a  typical  pattern  when  crossing  a  resonance  through  an  island. 


Figure  18:  Plots  of  rotation  number  for  different  cases.  Left  plot:  the  rotation  number  p  as  a  function  of 
the  initial  value  of  pr  with  r  =  r*  for  M  =  350  g  and  also  for  the  pulley-less  case  with  p  =  3.  A  zoom-in 
around  a  resonance  is  shown  in  the  inset.  Right  plot:  p  as  a  function  of  pr  for  M  =  180  g.  See  the  text 
for  additional  details. 


On  the  right  hand  side  of  Figure  [TH  we  show  the  results  for  M  =  180  g.  As  before, 
we  use  initial  points  (r*  ;  pr)  with  pr  going  away  from  the  fixed  point  of  Vm.  V ^  has 
also  been  used  instead  of  Tm  because  the  latter  is  close  to  a  symmetry  with  respect  to 
(r*  ;  p*),  hence  the  displayed  rotation  number  is  small.  It  increases  monotonously  as  we 
move  away  from  the  fixed  point,  also  marked  as  a  large  dot.  Now,  however,  the  intervals 
with  p  £  Q  constant  are  extremely  narrow.  We  observed  a  few  resonances,  such  as  2/45 
and  5/116,  checking  that  the  width  of  pr-intervals  is  below  10~6  kg  •  m  ■  s_1. 

We  want  to  mention  now  that  the  pattern  of  p  as  a  function  of  pr  is  a  clear  indication  of 
non-integrability.  Indeed,  let  us  first  look  at  the  left  plot  in  FigureHH  We  have  also  shown 
the  results  for  SAM  without  pulleys  and  p  =  3,  the  integrable  case.  The  level  of  energy  is 
the  one  corresponding  to  (ro  ;  9q)  =  (0.25  m ;  0.0°),  similar  to  the  kind  of  reference  orbits 
used  in  [3].  A  periodic  orbit  is  found  near  (r*  ;  p*)  =  (0.125  m,  0.369587  kg  •  m  •  s-1). 
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Using  r  =  r*  and  pr  between  p*  and  zero,  we  check  that  p  decreases  to  zero.  Under  V„, 
there  is  a  “separatrix”  bounded  by  r  =  0,  pr  =  0  and  a  curve  of  the  form  r  =  0.25  —  ap% 
for  a  suitable  a.  The  system  being  integrable,  the  map  V ^  has  a  first  integral  I  and  the 
iterates  are  on  level  curves  of  I.  When  approaching  the  separatrix,  the  dynamics  slows 
down  near  r  =  0  and  pr ,  hence  p  is  very  small. 

But  in  the  pulley  case,  the  rotation  number  ceases  to  exist  at  a  value  which  is  unrelated 
to  any  separatrix.  This  is  against  the  typical  behaviour  of  integrable  maps. 

One  should  expect  chaotic  dynamic  regions  in  a  resonance  zone  in  SAM  with  pulleys. 
Beyond  the  islands  around  periodic  elliptic  points,  there  are  homoclinic  tangles  associated 
to  the  hyperbolic  zones.  Attempts  to  visualize  them  lead  to  the  conclusion  that  the  size 
of  the  “chaotic  zones”  is,  at  most,  of  the  order  of  magnitude  of  the  round-off  errors  with 
standard  double  precision  computation.  Hence,  the  escape  remains  a  main  evidence  of 
non-integrability. 

6-4-  Remarks  on  dissipation 

As  we  have  already  noted,  experiments  show  non-negligible  dissipative  phenomena 
which  decrease  mechanical  energy  during  the  motion.  Hence,  equations  of  motion  such  as 
ED  which  do  not  contain  any  dissipation  term  yield  by  no  means  a  complete  description 
of  SAM  dynamics. 

In  order  to  show  that  the  observed  convergence  of  the  pendulum  towards  the  pulley 
(Figure  [7Jr)  is  due  to  dissipation  and  not  to  a  too  short  integration  time  interval,  the 
corresponding  theoretical  trajectory  (Figure  fl^hl  has  been  integrated  for  much  larger 
intervals.  Trajectories  are  not  shown  but  we  can  say  that  the  limited  region  of  space 
occupied  by  the  pendulum  in  Figure  I12h  is  progressively  filled  as  time  is  running.  The 
corresponding  Tufillaro  trajectory  integrated  for  larger  time  is  in  addition  symmetrical. 

Four  main  different  sources  of  dissipation  can  be  noted:  air  friction  on  the  pendulum 
and  the  counterweight,  thread  friction  on  the  pulley  grooves  due  to  an  imperfect  absence 
of  slippage  condition,  and  dissipation  due  to  the  ball  bearing  of  the  pulleys.  This  last 
source  is  implicitly  included  in  the  equation  of  motion  through  the  measure  of  Ip  since 
it  has  been  determined  using  a  pulley  (see  Sectional). 

An  estimation  of  air  friction  can  be  performed  by  comparing  the  weight  of  the  pen¬ 
dulum,  mg ,  to  the  drag  force  exerted  by  air  on  the  pendulum,  whose  expression  is: 

Fd  =  Cd*fs 

where  Cd  is  the  drag  factor,  pa  =  1.29  kg  •  m~3  the  air  density  at  ambient  temperature 
T  =  298  K,  v  the  velocity  of  the  ball,  and  S  the  effective  surface  of  the  pendulum.  For 
a  spherical  ball,  Cd  ss  0.4.  Since  v  ss  lm-s_1,  S  =  ttD^/A  with  Db  =  30mm,  the 
drag-weight  ratio  is  of  the  order  of  2  x  10-4.  One  gets  the  same  order  of  magnitude  for 
the  cylindrical  counterweight  for  which  Cd  ~  1.  Dissipation  will  be  the  detailed  topic  of 
a  further  paper. 
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7.  Non-integrability  of  SAM  with  pulleys 

In  this  Section,  a  rigorous  and  original  analytical  proof  of  non-integrability  of  SAM  is 
performed  in  order  to  complete  the  above  rotation  number  and  Poincare  section  analysis. 
Needless  to  say,  this  non-integrability  result  is  fundamentally  different,  both  in  approach 
and  scope,  from  the  numerical  and  graphical  inference. 

Since  SAM  may  be  expressed  in  terms  of  symplectic  formalism,  in  order  to  detect  or 
predict  chaotic  behaviour,  it  is  pertinent  to  recall  some  concepts  related  to  the  integra- 
bility  of  Hamiltonian  systems  in  the  sense  of  Liouville- Arnold. 

7.1.  Algebraic  background  for  studying  integrability 

7.1.1.  Linear  and  Hamiltonian  integrability 

Differential  Galois  Theory..  See  [Tsjj .  [32]  and  0  for  more  information.  Given  a  linear 
differential  system,  with  coefficients  in  a  differential  field  (. K ,  d)  whose  field  of  constants 
C  is  algebraically  closed  (e.g.  [C(f) ,  d/dt]), 

dy  =  A(t)y ,  (8) 

an  algebraic  group  G  exists,  called  the  differential  Galois  group  of  (0),  acting  over  the 
C- vector  space  (if1: . . . ,  xf>n)  of  solutions  of  (0)  as  a  linear  transformation  group  over  C. 
Furthermore,  G  contains  the  monodromy  group  of  ([5]).  The  Galoisian  formalism  proves 
useful  here  due  to  the  following:  (0)  is  integrable  (i.e.  its  general  solution  can  be  written 
as  a  finite  sequence  of  quadratures,  exponentials,  and  algebraic  functions)  if,  and  only  if, 
the  identity  component  G°  of  the  differential  Galois  group  G  of  0)  is  solvable. 

Everything  said  in  the  previous  paragraph  may  be  obtained,  mutatis  mutandis ,  for 
linear  homogeneous  differential  equations 

d™  dn_1  d 

an  (t)  -^y  +  a„— 1  (t)  d - \-ai(t)—y  +  a0  (f)  y  =  0. 

In  Subsection  17.1.31  we  will  denote  Galois  groups  in  this  Gal  (L)  accordingly,  L  £ 
C(t)  [d/dt]  being  the  corresponding  differential  operator.  See  lf§j  for  more  details. 

Liouville- Arnold  integrability..  On  the  other  hand,  we  call  a  Hamiltonian  system  q  = 
dH/dp ,  p  =  — dH/dq ,  whether  or  not  linear,  meromorphically  integrable  (in  the  sense  of 
Liouville- Arnold) ,  if  it  has  as  many  independent  meromorphic  first  integrals  in  pairwise 
involution  as  degrees  of  freedom.  Same  applies  in  the  above  definition,  mutatis  mutandis, 
substituting  algebraic,  rational  or  any  other  function  class  in  for  meromorphic.  For 
the  sake  of  simplicity,  conjugate  canonical  variables  will  be  henceforth  assembled  in  a 
single  vector  z  =  (q,p)  and  the  Hamiltonian  system  will  be  written  in  compact  form 
z  =  Xh(z).  Everything  is  considered  in  the  complex  analytical  setting  from  this  point 
on. 

7.1.2.  Morales- Ramis- Ziglin  Theory 

For  each  integral  curve  F  =  :  t  £  /|  of  a  given  autonomous  dynamic  system  in 

dimension  m 


z=X(z), 


(9) 
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the  variational  equations  of  order  k  for  (j9j)  along  F,  VEp,  are  satisfied  by  2^  :=  dk(j>/dzk 
~  see,  e.g.,  E3  for  explicit  expressions  of  VEp  for  general  k  in  terms  of  vectors  k  = 
(hi, . . . ,  km)  G  Z™  such  that  k  =  k\  +  . . .  +  km.  Note  that  dkcj)/dzk  is  an  abridged 
notation  for  the  (k  +  l)-dimensional  matrix  of  all  partial  derivatives  of  <f>.  In  other 
words,  it’s  a  vector  for  k  =  0,  a  matrix  for  k  =  1,  etc. 

In  particular,  for  k  =  1  and  denoting  the  matrix  of  the  first-order  variational  equations 
simply  by  2,  we  obtain 

3  =  X'  (0)  2.  (VEr) 

We  thus  have: 


•  a  (generally  nonlinear)  system  @  and 

•  a  linear  system  flVEpD  linked  to  ©. 


The  hallmark  theorem  in  this  approach  connects  the  two  notions  of  solvability  listed  in 
17.1.11  namely  as  applied  to  a  Hamiltonian  X 'h  and  its  variational  equations  VEp,  along 
an  integral  curve  F  of  Xh ■  The  whole  theory  is  actually  the  ad  hoc  implementation 
of  the  following  heuristic  principle:  if  a  Hamiltonian  is  integrable,  then  its  variational 
equations  must  also  be  integrable. 


We  assume  T,  a  Riemann  surface,  may  be  locally  parameterized  in  a  disc  /  C  C  of 
the  complex  plane;  we  may  now  complete  T  to  a  new  Riemann  surface  F,  as  detailed  in 
[3  §2.1]  (see  also  [3,  §2-3]),  by  adding  equilibrium  points,  singularities  of  the  vector 
field  and  possibly  t  =  00. 


Theorem  7.1  (J  Morales-Ruiz  &  J-P  Ramis,  2001).  Let  H  be  an  n- degree- of-freedom 
Hamiltonian  having  n  independent  first  integrals  in  pairwise  involution,  defined  on  a 
neighbourhood  of  an  integral  curve  T.  Then,  the  identity  component  Gal  (VEp)0  is  an 
abelian  group  ( i.e .  Gal  (VEp)  is  virtually  abelian/ 

See  [3;  Corollary  8]  or  [3,  Theorem  4.1]  for  a  precise  statement  and  a  proof. 


7.1.3.  Differential  operators.  A  primer  in  the  Boucher-Weil  Theorem 
Linear  differential  equations..  See  [3.  §2]  for  more  details.  Let 

)n  /  d  \  n—  1 

+  a-n- 1  (  J  +  •  •  ■  +  0,0  ;  an, . . . ,  ao  G  C  (t) 


be  a  differential  operator  with  coefficients  in  the  field  of  formal  Laurent  series.  If  L(y)  =  0 
has  a  solution  of  the  form  y  =  tp  Ylk>o  Cfc^  anc^  c°  ^  f°rmal  substitution  of  y 

into  the  differential  equation  yields  the  vanishing  of  all  powers  of  t,  the  smallest  one 
among  them  -  we  call  the  equation  P  (p)  =  0  derived  from  the  latter  vanishing  the 
indicial  equation  (at  0),  the  roots  of  which  are  usually  called  exponents  of  L  (y)  (at 
0).  In  particular  we  can  also  define  the  indicial  equation  at  infinity  by  means  of  the 
transformation  f  =  1/t  and  expansion  around  /  =  0. 

It  is  a  known  fact  ([3:  Lemma  2.1])  that  the  degree  of  P  (p)  =0  is  at  most  n.  A 
singular  point  c  of  L  is  called  regular  singular  if  deg  P  ( p )  =  n.  A  linear  differential 
equation  L  =  0  with  only  regular  singular  points  (including  00)  is  called  Fuchsian. 
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We  call  L  reducible  if  it  factors  in  a  product  of  operators  of  smaller  positive  order. 
An  operator  L  admits  a  first  order  factor  d/d t  —  /,  /  £  C  (£)  if  and  only  if  L  (y)  =0 
admits  a  solution  y  such  that  y  =  fy\  in  particular: 

Lemma  7.2.  If  L  is  of  order  2:  L  is  reducible  if,  and  only  if,  it  has  an  exponential 
solution,  i.e.  a  solution  whose  logarithmic  derivative  is  rational.  □ 

Lemma  7.3  (@,  §3.1.2]).  If  L  is  Fuchsian,  every  exponential  solution  must  be  of  the 
form  y  =  Ili=i  (t  ~  £j)6i P  (£) ,  where  P  £  C  [X]  and  t\ , . . . ,  ts  are  finite  singularities  of  L 
with  exponents  e±, . . . ,  es,  whether  or  not  integers. 

Normal  variational  equations.  Let 

T  _  (  0  Id„  \ 

^  — Idn  0  ) 

be  the  canonical  symplectic  matrix.  Given  a  Hamiltonian  system  z  =  Xh  (z)  =  JXH  ( z ) 
expressed  in  Darboux  canonical  coordinates  z  =  ( q ,  p) ,  system  flVErD  reads  5  =  JH"  (z)  5 
along  r  =  {z  (£)}. 

Gauge  transformations..  A  i— >  P  [A]  :=  P_1  (^AP  —  P^j,  P  being  a  given  symplectic 

matrix,  may  be  used  to  reduce  flVErD  by  selectively  vanishing  a  number  of  entries  in 
JH"  (z)  (see  §5.2],  [H,  §4.1],  p|)-  A  typical  first  choice  consists  of  symplectic  matrices 
of  the  form 

P  (.*)  =  (  C2  C3  C4),  (10) 

in  order  to  induce  a  row  and  a  column  of  zeroes  in  the  variational  matrix.  When  such  is 
the  case,  QVErP  acquires  a  “decoupled”  appearance  and  a  system  of  order  n—  2  may  be 
extracted  therefrom.  Such  a  system  is  usually  called  a  normal  variational  system  along 
r.  Let  us  denote  it  by  NVEr- 

As  is  always  the  case  with  all  differential  systems,  a  cyclic  vector  (0,  0)  may  be 
used  to  obtain  a  linear  differential  equation  L  (y)  =  0  of  order  n  —  2  equivalent  to  NVEp. 

The  result  central  to  this  paragraph,  and  a  particularly  useful  consequence  of  Theorem 
O  is  the  following  (see  also  0  Proposition  9  &  Theorem  8  (§5.3)]): 

Theorem  7.4  (D  Boucher  &;  J-A  Weil,  s  Criterion  1]).  Let  Xh  be  a  Hamiltonian 
system  and  L  its  normal  variational  operator  along  a  given  particular  solution.  If  L  is 
irreducible  and  displays  logarithms  in  a  formal  solution,  then  Gal(L)°  is  not  abelian,  i.e. 
Xh  is  not  integrable.  □ 

Typically,  the  most  difficult  part  in  trying  to  apply  Theorem  17.41  is  to  check  the 
irreducibility  condition. 


1.2.  Statement  of  the  main  result 

We  recall  Mt  =  M  +  m  +  2 Ip  /  R 2  and 


(p-2  +  Rpi)2 

mq\ 


+  gqi(M 


?7lCOS<72)  —  gR{Mq2 


msmq2), 


where  q\  =  r,  q2  =  9,  p\  =  pr ,  and  p2  =  pg.  The  main  result  in  this  Section  is  the 
following: 
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Theorem  7.5  (Non-integrability  of  SAM  with  massive  pulleys).  For  every  phys¬ 
ically  consistent  value  of  the  parameters,  regardless  of  Ip  and  R,  X-r  is  meromorphically 
non-integrable. 

This  is  a  complement  to  what  has  already  been  proved  for  SAM  without  pulleys ,  i.e. 
the  limit  case  Ip  =  0,  R  =  0  and  Mt  =  M  +  m: 

Uw  =  \  (lu  +  -^2)  +59i  iM  -  mcosq2) ;  (11) 

2  \Mt  mqf  J 

in  that  case,  the  following  held: 

Theorem  7.6  ( SAM  without  massive  pulleys). 

1.  (fS,  Theorem  1],  Equation  (16)],  0;  see  also  Remark [7T](2))  If  M  >  m  and 

_  M  p(p+  1) 

^  m  ^ p  p(p  +  1)  —  4 

for  every  p  £  h,p  >  2,  then  Hamiltonian  Xrw  is  non-integrable.  In  particular,  if 
p  £  (3/2,  3)  U  (3,  00),  it  is  non-integrable. 

2.  (0,  Equation  (16)])  For  p  =  2,  p  =  p2  =  3,  X-rw  is  integrable  and  has  the 
following  first  integral: 


I  =  Q1Q2  qic- 


qid2 


22  2  2  1  PldlC  2p2S 

+  gq1sc  =gq1cs+p2 — 4^2  - , 


where  c  =  cos  (52/2),  s  =  sin  ((72/2). 

3.  (0  Theorem  4])  The  degenerate  cases  pp,  p  >  2  referred  to  in  item  1  are  non- 
integrable.  □ 


The  last  case  is  significantly  more  difficult;  it  relies  on  the  higher  order  variational 
equations  17 1  and  uses  techniques  introduced  in  fl3|. 


7.3.  Proof  of  Theorem 

We  have  two  particular  solutions  for  SAM  with  massive  pulleys: 


zo,n(t)  = 


Ztt  .nid)  — 


g(M  -  m)  f  t(l  -t) 
2  \  Mt 

g(M  +  m)  / t(  1  —  t ) 


0,  1  -2t,  R(2t  -  1)  , 


2tt 


V  Mt  ’  g  (M  +  m) 


,1-2 1,  R(2t  -  1)  . 


System  (|VErD  around  V  =  {z^, /?(£)}  takes  the  form  5  =  Az.  with 

/0  0  1/Mt  +  R2ai  Ra  i\ 

0  0  Ra\  a\ 

A  = 


0  0  0 
V  0 


a  2 


0 


0 

0  / 
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where 

4  Mf  mg2(M  +  m)t(t  —  1) 

Ul  =  mg2(M  +  m)2t2(t  —  l)2  “  °2  =  2Mt  ' 

Notably,  flVErD  decouples  without  the  need  for  an  additional  gauge  transformation  such 
as  m  as  is  usual  and  customary  (0.0.  mi).  and  as  would  be  the  case  if  T  =  {zWio(t)}. 
See  Remark  1 7.  II 1  below.  System  NVEr  takes  the  form  =  B<J>  with 

bj°  “‘V 

V  «2  0  ) 


A  necessary  condition  of  integrability  is  the  virtual  abelianity  (see  Theorem  mi)  of 
Gal  (NVEr),  NVEr,  by  means  of  a  cyclic  vector  c  and  the  subsequent  gauge  transfor¬ 
mation  given  by  $  =  where  Q  =  (  c  c+BTc  )T  (jj,  §B.4]),  will  take  the 

form 


/ 

V 


o 

2  Mt 

( M  +  m)[t  —  1  )t 


1 

2(2t  —  1) 
(t  -  1  )t 


in  presence  of  a  constant  cyclic  vector  c  =  (a,  0).  The  above  system  is  obviously  equiva¬ 
lent  to  the  following  hypergeometric  (0,  §15.5]),  hence  Fuchsian  differential  operator: 


/d\2  2(2t  —  1)  d  2 Mt 

\d t)  ^  t(t  —  1)  dt  (M  +  m)t(t  —  1)  ’ 


Let  us  now  check  the  virtual  non-commutativity  of  Gal  ( L ).  L  has  three  singularities: 
0,  1,  oo.  At  t  =  0  or  t  =  1,  we  have  local  exponents  —1  and  0.  The  formal  solution  at 
t  =  0  has  logarithms  except  in  two  cases:  m  =  0,  Mt  =  M  +  m,  both  discarded  in  our 
case  since  they  would  correspond  to  no  small  mass  and  no  pulley,  respectively.  Indeed, 
a  particular  solution  is 


/  3  [9(M  +  m)  —  8-Mi]1/2  3  [9{M  +  m)  -  8Mt}1/2  .  n . 
\2  2(M  +  to)1/2  ’2+  2  (M  +  m)1/ 2  ’  ’ 


where 


2F1  ( a,b ;  c; 


OO 


*)  =  £ 

k= 0 


{a)k{b)k 
( c)k  k\ 


is  the  Gauss  hypergeometric  function  ( [T],  §15.1]),  ( a)k  =  a(a+l)  •  •  •  (a  +  fc—  1)  being  the 
Pochhammer  symbol.  This  first  solution  has  a  local  expansion  around  t  =  0  of  the  form 


2/1  =  1  + 


M  +  m 


t  + 


Mt  [2  (M  +  m)  +  Mt\  ^ 2  +  ^  \ 
3  (M  +  m)2  ^  1  ’ 


We  may  then  obtain  an  expansion  for  a  second  formal  solution 


1  Mt 
t  M  +  m 


2 Mt  [Mt  -  3 (M  +  m)\ 
( M  +  to)2 


_  /  n\  M  +  m  —  Mt  _  . 
t  +  O  ( t 2)  +  2 - — — - yi  Inf. 

M  +  TO 
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Keeping  the  Boucher- Weil  Theorem  17.41  in  mind,  and  in  presence  of  the  logarithm  in  1/2, 
there  is  obstruction  to  integrability  if  L  is  irreducible. 

Let  us  assume  it  is  reducible.  In  virtue  of  Lemma  O  L  admits  an  exponential 
solution.  We  recall  that  L  is  Fuchsian.  The  expansion  of  an  exponential  solution  around 
t  =  0  does  not  contain  logarithms,  although  we  have  shown  there  are  formal  solutions  with 
logarithms  around  the  singularities  -  as  has  been  seen  explicitly  for  t  =  0.  Thus,  those 
without  logarithms  correspond  to  the  maximal  exponents  ([3J]),  hence  the  admissible 
exponents  at  the  finite  singularities  are  all  0.  Hence,  in  virtue  of  Lemma  17.31  the  only 
possible  form  for  an  exponential  solution  y  would  be  that  of  a  polynomial  solution;  let 
N  be  its  degree.  Expanding  y  in  increasing  powers  of  i-1,  —  N  is  the  exponent  of  the 
leading  term,  hence  an  exponent  at  infinity:  oaktk  =  (l/i)_w  aN-k(^/t)k ■ 

Now,  the  exponents  at  infinity  are  the  two  roots  of  the  indicial  equation  p2  —  3p  + 
2 Mt/(M  +  to).  Since  —N  is  such  a  root,  this  means  Mt  =  —N(N  +  3 )(M  +  to)/ 2  with 
N  positive.  Therefore  Mt  would  be  negative,  which  is  physically  irrelevant. 

Hence  follows  that  the  Swinging  Atwood  Machine  system  with  massive  pulleys  is 
always  non-integrable  with  meromorphic  first  integrals.  □ 

Remarks  7.1. 

1.  Intriguingly,  the  solution  used  for  our  proof  was  which,  at  least  for  the  pulley¬ 
less  case  R  — >  0,  Mt  — >  M  +  to,  and  although  mathematically  plausible,  has  no 
actual  physical  significance.  The  solution  which  would  be  physically  possible  for  all 
values  of  R,  Zo„R,  posed  further  problems  with  regards  to  system  JVErD  and  was 
finally  discarded  in  our  proof.  It  is  worth  mentioning,  however,  that  in  the  case 
without  pulleys,  [l4j  used  precisely  the  corresponding  form  of  the  latter  solution 
zo,o- 

2.  The  same  proof  given  for  Theorem  17.51  may  be  obtained,  analogously,  for  SAM 
without  pulleys  T~LW ,  p  >  1 ,  and  the  proof  yields  precisely  item  1  in  Theorem  17. (>l 
Indeed,  by  using  the  “classical”  solution  z0, o  (which  corresponds  to  the  original 
Atwood  machine)  and  an  adequate  gauge  transformation,  we  obtain  the  normal 
variational  equation: 


2 Mt  .  ,  (2t  —  1)  d  d2  .  „ 

(M  -  m)t(t  - 1) y(t)  ~  (TVTjf  dAt)  +  de9*1*  =  °- 


This  is  a  Gauss  hypergeometric  equation.  Local  solutions  are 


f2- 


Mt 


3  V  M  —  to 


+  lj  t3  +  0  ( t 4) , 


and 


(M  -  to)  [2 t.Mt  +  {M  -  to)]  /  2n  , 
2 Mt  [Mt  +  (M  -  to)]  [  1 


t2- 


M, 


3  \  M  - 


+  1  )t3+o  (f4) 


Ini. 


The  degenerate  cases  Mt  =  0  and  Mt  +  M  —  m  =  0  are  of  course  not  physically 
acceptable.  The  exponents  at  zero  are  0  and  2,  the  exponents  at  t  =  1  are  also  0 
and  2,  and  the  exponents  at  infinity  are  the  roots  of  the  polynomial  Pco(X)  where 
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P00(X)  =  ( M  —  to)  X2  +  3  ( M  —  m)  X  —  2  Mt.  Same  as  in  the  proof  of  Theorem 
17.51  reducibility  would  occur  only  for  a  polynomial  solution,  i.e.  in  presence  of  an 
integer  N  such  that  Poo(N)  =  0,  implying 

Mt  =  J  (IV  +  4)  (IV  +  1)  (M  -  to)  . 

Setting  Mt  =  M  +  m  (that  is,  the  pulley-less  case)  we  would  have  an  equation  for 
p,  whose  solution  would  be 

N2  +  5N  +  6  (N  +  3)(N  +  2) 

M  "  N2  +  57V  +  2  "  (N  +  3)(N  +  2)  -  4’ 

obviously  equivalent  to  the  condition  in  item  1  in  Theorem  17.61  for  N  =  p  —  2. 

3.  For  f/,=  l,  however,  the  closest  thing  to  such  a  proof  as  that  sketched  in  item  2 
is  discarding  the  existence  of  first  integrals  with  a  meromorphic  growth  at  infinity, 
e.g.  rational  ones.  This  is  due  to  the  fact  that  the  normal  variational  equation 
around  particular  solution  z  (t)  =  (—gt/2  ;  0 ;  —  gm  ;  0)  is  a  Hamburger  equation , 


2  Id  d2 

=  °. 


i.e.  an  equation  with  exactly  two  singularities:  a  regular  one  at  zero  and  an  irregular 
one  at  infinity  [ll],  §17.6];  the  solutions  are  almost  Bessel  functions:  the  general 
solution  is: 

y(t)  =  Ci  1 1-2  (2y/2t)  +  c2t  K2  (2y/2t)  ,  Cl  ,  C2  g  C, 


where,  n  £  Z  given,  In  (2)  and  I\n  (z)  are  the  modified  Bessel  functions  of  the  first 
and  second  kind ,  respectively  (0P-  416],  0,  p.  376],  H,  p.  185]): 


1  f  exp  [(z/ 2)  (t  +  1  /t)\ 
2niJSi 


_  7T  I_n  (z)  -  In  jz) 
2  sin(ri7r) 


both  having  a  branch  cut  discontinuity  in  the  complex  z  plane  running  from  —00  to 
0  (although  In  is  regular  at  0,  whereas  Kn  has  a  logarithmic  divergence  at  0).  This 
assures  the  presence  of  a  non-trivial  Stokes  multiplier  at  infinity  for  the  variational 
equation  when  fi  =  1.  This  implies  the  following  important  conclusion:  the  system 
is  not  integrable  with  first  integrals  which  are  rational  functions  ofr,  9 ,  cos  9,  sin#, 
pr,  pg,  where  9  belongs  to  a  neighbourhood  of  0. 

Indeed,  any  first  integral  is  a  function  of  ( q,p ),  hence  must  be  defined  on  the 
phase  manifold  M  :=  C4.  We  can  partially  compactify  M  into  M  :=  P^.  x  C2, 
where  P^.  stands  for  the  compactification,  by  means  of  the  inverse  stereographic 
projection  P-1,  of  the  domain  C2  for  (qi ;  pi)  =  (r  ;  pr),  whereas  the  second  factor 
C2  is  the  (6  ;  pg  (-plane.  Using  P_1 in  order  to  compactify  the  particular  solution 
T  =  {  z  (t)}.  we  obtain  a  Riemann  sphere  I  in  P^.,  whose  immersion  in  P2  x  C2 
is  contained  in  {9  =  0};  therefore  cos#,  sin#  are  holomorphic,  hence  meromorphic 
on  a  neighbourhood  of  T  in  M.  Theorem  17.11  implies  the  absence  of  a  complete 
set  of  first  integrals  which  are  meromorphic  on  a  neighbourhood  of  T  in  the  partial 
compactification  M.  Since  any  rational  function  of  r,  9,pr,pg,  cos#,  sin#  must  be 
meromorphic  in  a  neighbourhood  of  T  in  M .  we  obtain  the  last  claim  in  the  previous 
paragraph. 
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8.  Conclusion  and  perspectives 

In  this  paper,  experimental  and  theoretical  results  concerning  the  Swinging  Atwood 
Machine  are  presented.  Equations  of  motion  with  two  pulleys  are  found  generalizing 
those  of  previous  studies.  An  experimental  device  for  SAM  is  constructed  and  described 
in  detail.  Experiments  are  conducted  and  the  trajectories  retrieved  from  computer  video 
analysis  closely  match  those  arising  from  numerical  solution  of  the  equations  of  motion. 
Such  comparisons  seem  to  show  that  the  motion  is  influenced  by  the  non-negligible 
dimension  and  the  rotation  of  the  pulleys;  in  particular,  a  non-zero  pulley  radius  leads  to 
asymmetric  pendulum  trajectories.  We  conclude  that  pulleys  cannot  be  ignored  in  SAM 
dynamics. 

Finally,  after  giving  some  numerical  evidence  of  the  lack  of  integrability  of  SAM  in 
the  absence  of  dissipation,  a  complete  proof  of  this  fact  is  given  using  differential  Galois 
theory  and  the  necessary  conditions  following  from  Morales-Ramis  theory. 

Several  perspectives  of  this  work  can  be  considered.  First  of  all,  other  experiments  are 
currently  conducted  with  /ieXp  =  2.03  and  qualitative  preliminary  results  (not  shown  in 
this  paper)  indicate  that  the  dynamics  of  SAM  seems  to  be  irregular  -  chaotic  behaviour 
is  expected.  Detailed  research  will  be  performed  in  a  future  work.  Secondly,  the  influence 
of  dissipative  phenomena  on  SAM  dynamics  should  be  studied;  a  possible  procedure  in 
such  direction  is  increasing  air  friction  by  coating  the  pendulum  with  different  materials, 
judiciously  chosen  so  as  to  induce  changes  in  drag  force.  Another  method  could  consist 
in  placing  the  counterweight  in  media  more  viscous  than  air,  such  as  water  or  glycerin, 
and  forcing  it  to  evolve  therein.  Perhaps  an  ultra  fast  camera  of  about  1000  images  per 
second  could  be  necessary  to  pick  up  much  more  points  and  achieve  a  better  resolution 
in  pendulum  trajectories. 

Concerning  the  integrability  of  SAM  without  pulleys,  that  is  Hamiltonian  Ttw  in  (fill) 
obtained  from  Q  by  skipping  the  contributions  of  R  and  Jp,  the  results  summarized  in 
Theorem  17.01  along  with  the  result  shown  in  the  previous  section  for  //  =  1 ,  close  the 
problem:  the  pulley-less  case  is  non-integrable  for  all  values  of  /r  7^  3. 
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